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SUMMARY 

The  Interactive  Zone  Embedding  concept  is  a  generalization  of  concepts 

which  have  been  used  for  some  time  with  boundary  layer  analyses.  This 

concept  has  been  extended  to  more  complex  flows  using  a  potential  flow  for 

the  outer  cone  and  the  compressible  Navier-Stokes  equations  for  the  inner 

zone.  In  this  interactive  manner  the  flow  about  a  rounded  trailing  edge 

was  computed  for  two  locations  of  the  outer  boundary.  The  calculations  were 

performed  on  a  5:1  elliptic  cylinder  at  a  Mach  number  of  0.2  and  a  Reynolds 

3 

number  based  on  chord  of  10  .  A  body  fitted  conformal  coordinate  system  was 
used.  In  the  first  instance  the  outer  boundary  was  located  approximately 
one  chord  from  the  ellipse  and  in  the  second  the  outer  boundary  was  brought 
to  about  one-half  chord  of  the  ellipse. 

The  displacement  body  was  chosen  as  the  interactive  streamline  and  the 
potential  flow  over  this  streamline  was  computed  as  the  sum  of  the  analytic 
potential  flow  over  the  elliptic  cylinder  plus  the  flow  from  a  series  of 
sources.  The  strengths  of  these  sources  were  determined  to  make  the  inter¬ 
action  streamline  also  be  a  streamline  in  the  potential  flowfield.  The 
resulting  potential  flow  was  used  to  evaluate  the  static  pressure  at  the 
outflow  boundary  of  the  Navier-Stokes  zone. 

The  calculations  were  performed  with  the  large  Navier-Stokes  domain 
and  were  repeated  with  a  smaller  Navier-Stokes  domain.  The  flow  properties 
in  the  vicinity  of  the  body  for  the  two  calculations  were  in  excellent 
agreement.  Farther  from  the  body  the  fluid  properties  were  in  somewhat 
poorer  agreement.  It  is  anticipated  that  refinements  in  the  interaction 
model  could  Improve  the  flowfield  prediction  away  from  the  body.  In  partic¬ 
ular,  it  is  felt  that  locating  the  inner  boundary  for  the  potential  flow 
calculation  further  from  the  body  and  taking  care  not  to  bring  the  downstream 
boundary  of  the  viscous  region  too  close  to  the  ellipse  trailing  edge  would 
be  very  beneficial. 


1 


INTRODUCTION 


The  cornerstone  of  modern  aerodynamics  is  the  ability  to  predict  the 
lift  and  drag  on  lifting  surfaces  such  as  airfoils,  or  cascades  of  airfoils 
when  concerned  with  turbomachinery  applications.  Early  aerodynamlcists 
recognized  that  at  high  Reynolds  numbers,  typical  of  aircraft  flight, 
viscous  effects  would  be  small  in  the  main  except,  as  Prandtl  pointed  out, 
very  close  to  the  airfoil  surface  within  the  boundary  layer.  With  this 
concept  it  was  possible  to  utilize  the  highly  developed  mathematical  theory 
of  potential  flows  to  predict  the  inviscid  irrotational  pressure  field 
acting  on  the  airfoil  or  cascade.  One  flaw  remained,  however,  and  that  was 
that  an  Infinity  of  solutions  could  be  obtained  from  potential  flow  theory 
and  that  some  physical  condition  had  to  be  invoked  to  determine  which 
solution  was  obtained  in  reality.  A  highly  satisfactory  resolution  to  this 
problem  for  incompressible  steady  flow  around  airfoils  with  sharp  trailing 
edges  was  developed  and  termed  the  Kutta-Joukowski  condition.  This  technique 
amounted  to  adjusting  the  circulation  about  the  airfoil  or  cascade  until  the 
rear  stagnation  point  occurred  precisely  at  the  sharp  trailing  edge.  Clearly 
such  a  technique  does  not  extend  to  rounded  trailing  edges  and  could  reason¬ 
ably  be  questioned  even  for  sharp  trailing  edges  in  unsteady  flow.  Given 
the  necessity  of  treating  potential  flows  through  cascades  of  airfoils  with 
rounded  trailing  edges,  the  turbomachinery  aerodynamlcists  developed  alter¬ 
native  strategies,  which  will  subsequently  be  discussed.  A  turbine  airfoil 
typical  of  modern  design  is  shown  in  Fig.  1. 

At  this  point  it  is  appropriate  to  digress  momentarily  and  discuss  the 
origin  of  rounded  trailing  edges.  In  external  or  wing  aerodynamics,  it  has 
usually  been  found  structurally  feasible  and  aerodynamlcally  desirable  (from 
a  drag  point  of  view)  to  maintain  a  sharp  trailing  edge  on  the  airfoil 
wherever  possible.  Some  advantage  may  be  gained  from  blunt  trailing  edges 
at  transonic  or  hypersonic  speeds  but  this  point  will  not  be  dwelt  on  here. 

In  turbomschlnes  blunt  trailing  edges  arise  frequently  simply  because 
menufacturlng  tolerances  do  not  permit  aerodynsmlcally  sharp  trailing  edges 
to  be  obtained.  Chord  lengths  of  one  inch  or  less  are  frequently  encountered 
in  turbomachines,  for  instsnce.  Secondly,  high  temperature  and  rotational 
environments  such  as  are  encountered  In  turbines  do  not  permit  the  use  of 


sharp  trailing  edges  due  to  resulting  stress  concentrations  which  would 
arise  at  these  sharp  points.  Trailing  edges  could  be  squared  off;  however. 
It  is  often  claimed  that  squared  off  trailing  edges  result  in  a  drag  penalty 
relative  to  rounded  trailing  edges,  and  consequently  rounded  trailing  edges 
have  retained  their  popularity  over  the  years.  Typically  trailing  edge 
radii  can  be  about  one  quarter  of  the  leading  edge  radius  and  in  some 
applications,  particularly  small  gas  turbines,  may  even  approach  the  leading 
edge  radius  In  size. 

Returning  now  to  the  problem  of  determining  a  suitable  alternative  to 
the  Kutta-Joukowski  condition  for  rounded  trailing  edges,  the  procedure 
presently  in  widespread  use  for  steady  flow  is  based  upon  cascade  data 
correlations.  The  overall  approach  and  one  data  correlation  is  given  in 
detail  by  Lleblein  (Ref.  1).  The  essence  of  the  procedure  is  that  the 
Kutta-Joukowski  condition  is  replaced  by  a  specified  turning  angle  of  the 
flow  through  the  cascade  obtained  from  a  data  correlation.  Camber,  solidity 
and  air  inlet  angle  being  the  essential  parameters  of  the  correlation.  The 
correlation  is  usually  presented  in  terms  of  a  deviation  angle  which  is 
defined  as  the  difference  between  the  average  direction  of  the  exiting  flow 
and  the  direction  of  the  blade  mean  line  at  the  trailing  edge.  Lleblein 
(Ref.  1)  cites  the  work  of  Carter  and  Hughes  (Ref.  2)  as  providing  a  theo¬ 
retical  basis  for  such  a  correlation,  and  much  subsequent  work  by  industry 
has  concentrated  on  expanding  the  correlation  data  base,  at  considerable 
cost  and  with  the  obvious  problems  of  extrapolation  to  new  designs  and 
different  flow  regimes.  The  steady  and  unsteady  transonic  cases  are,  to  say 
the  least,  somewhat  clouded  at  the  present  time.  Recently,  Gostelow  (Ref.  3) 
has  reviewed  the  question  of  the  location  of  the  trailing  edge  stagnation 
point  for  rounded  trailing  edges  and  suggested  a  rather  Intuitive  extra- 


1.  Lleblein,  S.:  Experimental  Flow  in  Two-Dimensional  Cascades.  Chapter 
VI,  NASA  SP  36,  Ed.,  I.  A.  Johnsen  and  R.  0.  Bullock,  1965. 

2.  Carter,  A.  D.  S.  and  Hughes,  H.  P.s  A  Theoretical  Investigation  Into  the 
Effect  of  Profile  Shape  on  the  Performance  of  Aerofoils  in  Cascade. 

Aero.  Research  Council  RAM  2384,  1946. 

3.  Gostelow,  J.  P.s  Trailing  Edge  Flows  Over  Turbomachine  Blades  and  the 
Kutta-Joukowski  Condition.  University  of  Cambridge,  Dept,  of  Engineering 
Report  CUED/ATurbo/TR  55,  1974. 


polacion  procedure  to  replace  the  deviation  angle  correlation  procedure 
currently  in  widespread  use.  Gostelow's  suggestion,  typical  of  the  present 
level  <*f  understanding  of  this  problem,  is  simply  to  perform  a  number  of 
potential  flow  calculations  with  assumed  deviation  angles  and  select  the 

deviation  angle  as  being  that  angle  which  allowed  the  upper  and  lower 
surface  trailing  edge  pressures  to  meet  when  the  pressure  is  linearly  extra¬ 
polated  from  the  85Z  chord  point.  Subsequently,  Miller  and  Serovy  (Ref.  4) 
reviewed  four  of  the  more  recent  suggested  methods  of  estimating  deviation 
angles  and  observed  that  the  suggestion  of  Gostelow  (Ref.  3)  and  the 
suggestion  of  Miller  and  Servoy  (Ref.  4)  both  gave  comparable  and  fairly 
good  results  for  a  very  limited  sample  of  data  near  the  design  Incidence. 

Off  design,  the  predictions  were  not  in  particularly  good  agreement  with 
the  measured  deviation  angles.  In  all  the  techniques  examined,  the  proposed 
methods  of  obtaining  the  deviation  angle  were  rather  intuitive  and  not  at 
all  rigorously  based.  Miller  and  Serovy  (Ref.  4)  concluded  that  "The  general 
lack  of  agreement  between  estimated  and  measured  deviation  angles  using 
inviscid  calculations  suggests  that  viscous  effects  must  be  Included  in  any 
general  deviation  angle  estimation  procedure."  Gostelow  (Ref.  3)  concluded 
that  for  the  three  suggested  deviation  angle  estimating  procedures  he 
examined,  including  his  own  and  the  Miller-Serovy  (Ref.  4)  suggestion,  that 
"All  three  procedures  fail  under  conditions  of  high  Incidence  or  loading." 

The  published  discussion  following  the  Miller-Serovy  paper  (Ref.  4)  clearly 
indicates  a  consensus  of  opinion  favoring  the  development  of  a  more  rigorous 
deviation  angle  prediction  scheme  which  would  Incorporate  viscous  and 
eventually  transonic  effects,  thus  providing  much  of  the  impetus  behind  the 
present  proposal.  Earlier  attempts  to  Include  viscous  effects  such  as  that 
due  to  Gostelow  et  al.,  (Ref.  5)  were  restricted  to  unseparated  trailing 


4.  Miller,  M.  J.  and  Serovy,  G.  K. :  Deviation  Angle  Estimation  for  Axial 
Flow  Compressors  Using  Inviscid  Flow  Solutions.  Trans,  of  ASME  J.  of 
Engineering  for  Power,  April  1975. 

5.  Gostelow,  J.  P.,  Lewkowicz,  A.  K.  and  Shaalan,  M.  R. :  Viscosity  Effects 
on  the  Two-Dimensional  Flow  in  Cascades.  Aeronautical  Research  Council 
of  Great  Britain  Current  Paper  CP  872,  1967. 


edges  or  used  highly  simplified  treatments  of  the  separated  flow  region 
(Geller,  Ref.  6). 

In  order  to  demonstrate  the  importance  of  the  stagnation  point  location 
(or  equivalently,  the  deviation  angle)  Gostelow  (Ref.  3)  presents  three 
calculated  pressure  distributions  around  a  cascade  of  airfoils  with  rounded 
trailing  edges  for  minor  changes  in  the  stagnation  point  location.  The 
resulting  changes  in  the  predicted  pressure  distribution  are  unfortunately 
very  significant.  Gostelow' s  results  are  reproduced  in  Fig.  2.  Also  taken 
from  Gostelow' s  report  are  the  typical  predicted  deviation  angles  which 
follow  from  the  various  closure  hypotheses  as  compared  with  measured  values 
for  a  particular  cascade  and  these  results  are  given  here  in  Fig.  3.  Mot 
mentioned  in  these  more  recent  studies  on  deviation  angle  prediction  is  the 
concurrent  very  poor  situation  which  must  exist  with  regard  to  loss  predic¬ 
tions.  One  bright  aspect  is,  as  Lleblein  (Ref.  1)  points  out,  that  the 
effect  of  compressibility  on  the  deviation  angle  correlations  is  not  large 
at  least  to  incidences  below  the  first  appearance  of  sonic  velocity. 

Clearly,  however,  transonic  flow  must  have  a  large  and  poorly  understood 
effect  on  the  flow  in  the  region  of  the  rounded  trailing  edge. 

Before  leaving  this  section  and  with  the  foregoing  background  it  is 
probably  appropriate  now  to  define  the  terminology  "rounded  trailing  edge" 
a  little  more  precisely  since  on  a  microscopic  level  aerodynamically  "sharp" 
trailing  edges  could  be  regarded  as  blunt.  Two  scales  enter  the  problem, 
the  first  an  invisdd  one  and  following  the  examples  shown  by  Gostelow 
(Ref.  3)  for  instance,  uncertainties  in  the  stagnation  point  location  of 
order  0.5Z  chord  create  significant  changes  in  the  predicted  pressure  field. 
Clearly  the  trailing  edge  radii  would  have  to  be  much  less  than  0.2SZ  chord 
to  tie  the  stagnation  point  down  sufficiently  well  that  the  airfoil  trailing 
edge  could  be  termed  "sharp".  As  a  working  hypothesis  in  the  subsequent 
work  a  trailing  edge  radius  of  less  than  0.02SZ  chord  will  be  considered 

6.  Geller,  W. :  Incompressible  Flow  Through  Cascades  With  Separation. 

AGARDograph  No.  164,  Boundary  Layer  Effects  in  Turbomachines,  Ed., 

J.  Surugue,  1972. 


"sharp”.  Such  a  value  could  be  expected  to  vary  somewhat  with  airfoil  shape 


but  serves  the  present  purposes  adequately.  Experimentally  Mash  et  al., 
(Ref.  7)  and  Inoue  and  Kaneko  (Ref.  8)  for  instance,  found  that  reducing 
the  trailing  edge  thickness  below  about  2Z  chord  had  little  effect  on  the 
observed  airfoil  or  cascade  performance,  respectively.  This  experimental 
observation  unfortunately  does  not  translate  into  the  fact  that  the  precise 
.  location  of  the  stagnation  streamline,  or  equivalently  the  deviation  angle, 
is  unimportant  in  calculation  schemes  when  the  trailing  edges  are  thinner 
than  2Z  chord.  Gostelow  (Ref.  3)  has  clearly  demonstrated  the  unfortunate 
sensitivity  of  the  flow  predictions  for  much  thinner  trailing  edges  than 
2Z  chord. 


7.  Nash,  J.  F. ,  Qulncey,  V.  C.,  and  Callinan,  J. :  Experiments  on  Two- 
Dimensional  Base  Flow  and  Subsonic  and  Transonic  Speeds.  Aero.  Research 
Council  of  Great  Britain  R&M  3427,  1966. 

8.  Inoue,  M.  and  Kaneko,  K. :  Effect  of  Trailing  Edge  Thickness  on  the 
Cascade  Performance  of  Circular-Arc  Blade.  2nd  Int.  JSME  Symposium 
Proceedings,  Tokyo,  September  1972. 


ANALYSIS 


Overview 

This  study  Includes  the  calculation  of  the  flow  In  the  vicinity  of  the 
founded  trailing  edge  of  an  airfoil.  In  performing  flow  field  calculations 
there  is  a  constant  desire  to  improve  efficiency  by  limiting  the  computa¬ 
tional  domain  such  that  the  available  grid  points  are  used  only  where  they 
are  absolutely  necessary.  Far  field  boundary  conditions  for  the  computa¬ 
tional  domain  are  rigorously  known  only  at  infinity.  Although  transformations 
are  available  to  map  grid  points  to  infinity,  greater  efficiency  can  be 
achieved  by  dividing  the  flow  field  into  inner  and  outer  zones  and  by  using 
computational  analyses  specially  suited  to  each  zone.  Prandtl's  boundary 
layer  theory  is  a  remarkably  successful  example  of  this  philosophy.  The 
outer  layer  or  zone  is  presumed  to  be  lnviscld  and  is  usually  treated  as  a 
potential  flow  which  is  readily  predicted  for  the  given  body.  The  boundary 
layer  or  inner  zone  is  then  computed  with  boundary  conditions  representing 
the  imposed  inviscid  pressure  gradients.  Although  in  some  cases  interaction 
between  zones  is  negligible,  there  are  many  instances,  especially  in  the 
cases  of  flows  containing  separation,  where  the  interaction  must  be  modeled 
to  maintain  accuracy.  The  general  technique  of  Interactive  Zone  Embedding 
is  directly  extendable  to  multi-zone,  three-dimensional,  and  time  dependent 
problems. 

Related  Work 

Three  decades  ago  a  great  deal  of  progress  was  made  in  understanding 
the  interaction  of  co-flowing  streams.  Tsien  and  Fins ton,  Ref.  9,  and 
Lighthlll,  Ref  10,  studied  the  interaction  of  an  inviscid  subsonic  stream 


9.  Tsien,  H.  S.  and  Finston,  M. :  Interaction  Between  Parallel  Streams  of 
Subsonic  and  Supersonic  Velocities.  Journal  of  the  Aeronautical 
Sciences,  Vol.  16,  No.  9,  Sept.,  1949. 

10.  Lighthlll,  M.  6. :  Reflection  at  a  Laminar  Boundary  Layer  of  a  Weak, 
Steady  Dlstrubance  to  a  Supersonic  Stream,  Neglecting  Viscosity  and 
Heat  Conduction.  Quarterly  Journal  of  Mechanics  and  Applied  Mathematics, 
Vol.  Ill,  1950. 
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along  a  surface  and  an  outer  supersonic  stream.  Crocco  and  Lees  (Ref.  11) 
extended  this  work  by  studying  the  Interaction  of  a  boundary  layer  and  an 
outer  supersonic  stream.  The  boundary  layer  was  treated  as  quasi-one- 
dimensional  using  Integral  boundary  layer  theory  and  the  supersonic  flow  was 
modeled  by  the  Prandtl-Meyer  relation.  Using  these  simple  analyses  for  the 
two  streams,  the  rate  of  growth  of  the  displacement  thickness  was  used  to 
determine  the  outer  flow  angle.  This  angle  In  turn  was  used  to  set  the 
boundary  layer  static  pressure.  This  technique  was  used  to  compute  flows  in 
adverse  pressure  gradients  up  to  the  point  of  separation. 

The  zone  interaction  concept  has  received  considerable  impetus  as  efforts 
have  been  made  to  extend  boundary  layer  theory  into  the  separated  flow  region. 
Studies  by  Stewartson  (Refs.  12,  13)  and  others  developed  triple  deck,  where 
the  boundary  layer  -is  divided  into  an  upper  layer  which  is  rotational  but 
inviscid  and  a  lower  layer  which  can  be  described  by  the  conventional 
incompressible  boundary  layer  equations.  The  outer  third  layer  or  "free 
stream"  is  perturbed  by  the  boundary  layer  presence  according  to  Ligh thill's 
concept. (Ref .  14).  The  three  layers  are  linked  by  the  pressure  field  and  an 
Iteration  devised  to  develop  the  mutually  compatible  interaction.  Much  of 
interest  about  separation  and  Interaction  can  be  learned  from  triple  deck, 
but  as  a  predictive  tool  it  is  restricted  to  asymtotically  high  Reynolds 
numbers  and  an  inner  layer  where  the  boundary  layer  approximations  are  valid. 
The  relationship  between  triple  deck  and  numerical  studies  of  interacting 
boundary  layers  has  been  examined  by  Davis  and  Werle  (Ref .  15) .  In  the 

11.  Crocco,  L.  and  Lees,  L.:  A  Mixing  Theory  for  the  Interaction  Between 
Dissipative  Flows  and  Nearly  Isentropic  Streams.  Journal  of  the 
Aeronautical  Sciences,  Vol.  19,  No.  10,  Oct.,  1952. 

12.  Stewartson,  K. :  Multlstructural  Boundary  Layers.  Advances  in  App. 

Mech.,  Vol.  14,  1974,  pp.  145-239. 

13.  Stewartson,  K. ,  Williams,  P.  G.:  Self  Induced  Separation.  Proc.  Roy. 
Soc.  Ser.A,  Vol.  312,  1969,  pp.  181-206. 

14.  Lighthill,  M.J.:  On  Boundary  Layers  and  Upstream  Influence.  Proc. 

Roy.  Soc.  Ser.  A,  Vol.  217,  1953,  pp.  478-507. 

15.  Davis,  R.  T.,  Werle,  M.  J. :  Numerical  Methods  for  Interacting  Boundary 
Layers.  Proceedings  of  the  1976  Heat  Transfer  and  Fluid  Mechanics 
Institute,  Stanford  University  Press,  1976,  pp.  317-339. 


conventional  interacting  boundary  layer  approach,  essentially  the  lower  and 
middle  decks  of  triple  deck  are  merged  and  solved  numerically,  together  with 
an  upper  deck  which  represents  the  interaction,  and  the  result  has  been 
termed  "composite  equations".  A  similar  approach  was  taken  by  Briley  and 
McDonald  (Ref.  16),  using  both  the  Navier-Stokes  and  boundary  layer  equations 
for  the  inner  region  and  a  potential  flow  interaction  model  to  provide  outer 
boundary  conditions.  By  comparison  of  interacted  Navier-Stokes  and  interacted 
boundary  layer  solutions,  it  was  shown  in  Ref.  16  that  for  small  thin 
separation  bubbles  the  conventional  boundary  layer  approximations  remain 
valid  when  streamwise  "elliptlcity"  arising  both  from  reverse  flow  velocities 
and  subsonic  displacement  Interaction  effects  are  accounted  for.  Two  general 
approaches  for  computing  zone  embedding  flow  models  with  interaction  have 
been  suggested  and  both  are  reviewed  from  the  viewpoint  of  boundary  layer 
separation  by  Williams  (Ref.  17).  In  the  direct  approach  (Ref.  16),  a  time- 
dependent  Iteration  path  is  followed  to  obtain  interacted  velocity  (or 
pressure  gradient)  boundary  conditions  for  the  inner  region  which  match  the 
outer  inviscid  flow  corrected  for  displacement  Interaction  with  the  inner 
viscous  layer.  The  inner  layer  equations  are  solved  as  a  time-dependent 
boundary  layer  or  Navier-Stokes  problem.  A  pressure  correction  to  some 
unperturbed  inviscid  body  pressure  distribution  is  computed  by  a  (transient) 
displacement  surface  Interaction  and  is  imposed  as  a  direct  boundary  condition 
for  the  viscous  layer.  The  computational  mesh  is  adjusted  at  each  time  step 
to  encompass  the  separation  region  both  upstream  and  downstream  and  includes 
the  transiently  changing  boundary  layer  thickness.  In  Ref.  16,  using  the 
boundary  layer  approximations,  downstream  boundary  conditions  or  vortlcity 
are  allowed  for  by  taking  implicit  sweeps  normal  to  the  wall  and  parallel 
to  the  wall  at  each  time  step.  The  composite  system  is  then  integrated  out 
to  steady  state. 

16.  Briley,  W.  R. ,  McDonald,  H.:  Numerical  Predictions  of  Incompressible 

Separation  Bubbles.  J.  Fluid  Mech.,  Vol.  69,  Part  4,  1975,  pp.  631-656. 

17.  Williams,  J.  C.,  Ill:  Incompressible  Boundary-Layer  Separation. 

Ann.  Rev.  Fluid  Mechanics,  1977. 
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In  the  interacted  inverse  approach  Carter  (Ref.  18)  assumed  a  distri¬ 
bution  of  some  quantity,  in  this  instance  displacement  thickness,  and  solved 
the  (steady)  boundary  layer  equations  for  the  imposed  pressure  gradient.  The 
difference  between  the  computed  pressure  gradients  obtained  from  inviscid 
flow  concepts  using  the  assumed  or  computed  boundary  layer  displacement 
distribution  and  the  pressure  gradients  determined  from  the  boundary  layer 
calculations  is  used  to  correct  the  assumed  skin  friction  or  displacement 
thickness  distribution.  Since  in  concept  the  direct  (Ref.  16)  and  indirect 
(Ref.  18)  approaches  are  merely  different  strategies  using  the  same  general 
displacement  surface  pressure  interaction  approach  to  obtain  solutions  for 
essentially  the  same  steady  equations,  there  should  in  principle  be  little 
difference  in  the  solutions  generated.  The  observation  can  be  made  however, 
that  this  Indirect  approach  may  Involve  spatial  marching  of  the  boundary 
layer  into  a  region  of  separated  flow  and  hence  either  an  (expensive) 
additional  Iteration  must  be  undertaken  to  allow  for  the  reverse  flow 
(Refs.  18,  19),  or  an  additional  approximation  Invoked  such  as  "FLARE", 
neither  of  which  are  required  in  the  direct  approach.  FLARE,  an  acronym 
for  Flugge-Lotz  and  Reyner  (Ref.  20),  stabilizes  the  computation  by  taking 
a  small  fraction  of  the  absolute  value  of  the  axial  convective  term  in  the 
reversed  flow  in  order  to  forward  march  into  the  separated  region.  As  might 
be  expected.  Carter  (Ref.  18)  found  that  as  the  recirculation  velocities 
decreased  the  differences  between  reverse  flow  iteration  and  FLARE  decreased 
for  thin  separation  bubbles  using  the  Interacted  inverse  boundary  layer 
approach.  As  a  corollary  to  the  similarity  of  the  steady  state  equations 
solved  by  the  direct  and  indirect  approaches  (FLARE  excluded),  it  follows 
that  results  obtained  with  the  Interacted  Inverse  approach  (Ref.  18)  which 
Indicate  that  it  deals  adequately  with  the  separation  singularity  problem 
should  be  equally  valid  for  the  Interacted  direct  approach  (Ref.  16). 

18.  Carter,  J.  E. :  Inverse  Solutions  for  Laminar  Boundary-Layer  Flows  with 
Separation  and  Reattachment.  NASA  Technical  Report,  NASA  TR  R-447, 
November,  1975. 

19.  Kllneberg,  J.  M. ,  Stager,  J.  L. :  On  Laminar  Boundary-Layer  Separation. 
AIAA  Paper  74-94,  1974. 

20.  Rayhner,  T.  A.,  Flugge-Lotz,  I.:  The  Interaction  of  a  Shock  Wave  With  a 
Laminar  Boundary  Layer.  Int.  J.  Non-Linear  Mech.,  Vol.  3,  No.  2, 

June  1968,  pp.  173-199. 
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In  the  case  of  supersonic  external  flow  Levy,  Shamroth,  Clbellng  and 
McDonald  (Ref.  21)  solved  a  shock  wave  turbulent  boundary  layer  Inter¬ 
action  problea  by  an  interaction  technique.  The  inner  zone  was  computed 
by  the  steady  finite  differenced  boundary  layer  equations.  The  small 
magnitude  of  the  reversed  flow  velocities  permitted  the  FLARE  stabilization 
to  be  used  without  excessive  error.  The  outer  zone  was  modeled  by  the 
Prandtl-Meyer  relation  applied  at  the  displacement  body  as  the  boundary  layer 
was  marched  downstream.  Since  neither  the  forward  marching  boundary  layer 
solution  nor  the  Prandtl-Meyer  relation  is  elliptic,  a  shooting  technique 
was  used  to  provide  the  upstream  influence  of  the  shock-boundary  layer 
Interaction.  The  calculations  of  Ref.  21  were  in  excellent  agreement  with 
the  experimental  data  of  Spaid  and  Frishett  (Ref.  22). 

Interactive  Zone  Embedding 

Zone  Embedding  is  a  technique  for  increasing  computational  efficiency 
by  dividing  the  flow  field  into  two  or  more  regions  and  using  appropriate 
numerical  analyses  for  each  region.  An  example  of  zone  embedding  is  the 
computation  of  flow  about  an  airfoil  by  solving  the  Navier-Stokes  equations 
in  the  vicinity  of  the  wing  with  the  farfield  boundary  conditions  determined 
by  potential  flow  about  the  airfoil. 

Interactive  Zone  Embedding  improves  the  accuracy  of  zone  embedding  by 
accounting  for  the  effect  the  fluid  dynamic  processes  in  each  zone  have  on 
the  fluid  in  adjacent  zones.  Three  steps  are  required  to  Implement  inter¬ 
active  zone  embedding.  The  first  step  is  to  choose  the  analyses  for  the 
different  parts  of  the  flow.  In  regions  where  the  stresses  are  important 
they  of  course  must  be  included  in  the  analysis.  However,  in  regions  where 
the  stress  is  negligible,  a  simplified  analysis  is  often  justified.  Knowing 
the  analyses  to  be  used,  the  second  step  is  the  division  of  the  flow  field 


21.  Levy,  R. ,  Shamroth,  S.  J. ,  Glbeling,  H.  J.,  and  McDonald,  H.:  A  Study 
of  the  Turbulent  Shock  Wave  Boundary  Layer  Interaction,  AFFDL-TR-7 6-163, 

-  Feb.  1977. 

22.  Spaid,  F.  W.,  and  Frishett,  J.  C. :  Incipient  Separation  of  a  Super¬ 
sonic  Turbulent  Boundary  Layer,  Including  Effects  of  Heat  Transfer. 

AIAA  Journal,  Vol.  10,  No.  7,  July  1972. 
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into  zones.  The  flow  in  each  of  these  zones  is  analyzable  by  a  different 
method.  The  determination  of  the  zones  is  subjective  and  care  must  be 
exercised  in  maintaining  consistency  with  the  fluid  dynamic  analyses  to 
be  used.  The  third  and  equally  crucial  step  is  the  determination  of  coupling 
relations  to  carry  Information  between  adjacent  zones.  Each  pair  of  adjacent 
zones  requires  a  pair  of  relations  to  communicate  the  effect  of  each  zone 
on  its  neighbor.  The  flow  of  information  provided  by  the  coupling  relations 
is  shown  schematically  in  Fig.  4.  Since  the  adjacent  analyses  do  not  all 
contain  the  same  amount  of  information  about  the  fluid,  the  coupling 
relations  can  transmit  only  the  "least  common  denominator"  of  Information 
contained  in  adjacent  analyses.  For  example,  modeling  the  displacement  body 
influence  on  external  flow  does  not  model  the  heat  transfer  or  the  shear  at 
the  Interface.  If  these  additional  properties  are  negligible  at  this  part 
of  the  flow  field  and  if  they  are  not  contained  in  the  external  flow 
analysis,  the  displacement  modeling  may  be  well  chosen.  If  these  properties 
are  modeled  in  the  external  flow,  their  values  at  the  interface  should  be 
included  in  the  interaction  model.  These  coupling  relations  contain  the 
exchange  of  information  -  momentum,  mass,  energy,  vorticlty,  etc.,  across 
an  Interface. 

Interactive  Zone  Embedding  is  directly  applicable  to  three-dimensional 
and  time-dependent  problems.  Navler-Stokes  calculations  about  finite  span 
wings  could  use  this  approach  to  limit  the  computational  domain  to  the 
vicinity  of  the  surface.  Oscillating  airfoil  and  cascade  calculations  could 
limit  the  Navier-Stokes  zone  to  the  vicinity  of  the  airfoils  by  using  a  time- 
dependent  potential  flow  analysis  to  model  the  far  field  wave  effects. 

Clearly  the  payoff  for  Interactive  Zone  Embedding  comes  from  using  a  very 
inexpensive  outer  analysis  to  restrict  the  relatively  expensive  analysis 
to  a  confined  region. 


APPLICATION  TO  FLOW  AT  A  ROUNDED  TRAILING  EDGE 


Interactive  Zone  Embedding,  as  presented  in  the  previous  section,  was 
used  to  analyze  the  flow  at  a  rounded  trailing  edge.  The  Navier-Stokes 
equations  were  used  in  the  zone  near  the  airfoil  since  the  boundary  layer 
equations  were  judged  to  be  inadequate  for  the  relatively  large  recirculation 
regions  expected.  In  the  outer  zone  a  potential  flow  was  used.  The  coupling 
relation  from  the  inner  zone  to  the  potential  flow  was  the  specification  of 
a  streamline,  and  the  coupling  relation  from  the  outer  zone  to  the  Navier- 
Stokes  was  the  specification  of  static  pressure.  The  actualization  of  this 
analysis  will  now  be  discussed. 

To  simplify  the  problem  during  the  initial  phase  of  the  study,  the  flow 
over  the  rear  of  a  5:1  elliptic  cylinder  was  considered.  An  orthogonal 
conformal  body  fitted  coordinate  system  was  used  consisting  of  ellipses  in 
the  azimuth  and  hyperbolas  in  the  radial  direction  as  shown  in  Fig.  5.  Since 
these  initial  calculations  were  run  at  zero  angle  of  attack,  syonetrlc  flow 
was  assumed  and  only  the  rear  upper  quadrant  of  the  flow  field  was  computed. 
The  equations  relating  the  computational  coordinates  y^,  y2  to  the  physical 
plane  x^,  x2  are 

x,  •  cosh  (y,)  cos(y2) 
x,  ■  sinh  (y.)  sin(y#) 

*  1  *  (l) 

Since  these  coordinates  are  conformal  thare  is  only  one  metric,  h,  at  each 
grid  point.  Sine  h  must  satisfy 


/ hdy(  ■  S, ,  /hdyt-  Sz 


(2) 


where  S^,  S2  are  the  physical  distances  along  the  coordinate  lines,  h  can  be 
determined  from  either 


h 


l/C 


as 

h  •  cosh*  Cy, I  -  cos*(yt) 


(3) 


where  and  are  Che  velocity  components  along  Che  coordinate  lines,  0  is 
Che  density,  p  is  the  static  pressure,  p,  p,  and  u  are  reference  pressure 
density  and  velocity.  Re  is  the  reference  Reynolds  number,  and  and  Fj  are 
the  viscous  force  components  of: 

F  ■?•["+  "*] 

where  7  and  7  t  are  the  viscous  stress  tensor  and  the  turbulent  average  stress 
tensor  respectively. 

These  equations  are  solved  by  the  linearised  Bock  Implicit  technique  of 
Briley  and  McDonald  (Ref.  23).  The  governing  equations  are  replaced  by  an 

23.  Briley,  W.  R.,  and  McDonald,  H. :  Solution  of  the  Multidimensional 

Compressible  Navler-Stokes  Equations  by  a  Generalised  Implicit  Method. 
Journal  of  Computational  Physics,  Vol.  24,  Aug.  1977,  p.  372. 


Implicit  time  difference  approximation,  optionally  a  backward  difference 
or  Crank-Nicoleon  acheae.  Terms  involving  nonlinearities  at  the  laplicit 
time  level  are  linearized  by  Taylor  expansion  about  the  solution  at  the  known 
tlae  level,  and  spatial  difference  approximations  are  introduced.  The 
result  is  a  system  of  multidimensional  coupled  (but  linear)  difference 
equations  for  the  dependent  variables  at  the  unknown  or  implicit  time  level. 

To  solve  these  difference  equations,  the  Douglas-Gunn  (Ref.  24)  procedure 
for  generating  alternating-direction  implicit  (ADI)  schemes  as  perturbations 
of  fundamental  implicit  difference  schemes  is  introduced  in  its  natural 
extension  to  systems  of  partial  differential  equations.  This  technique 
leads  to  systems  of  coupled  linear  difference  equations  having  narrow  block- 
banded  matrix  structures  which  are  solved  efficiently  by  standard  block- 
elimination  methods. 

Symmetry  boundary  conditions  were  imposed  on  the  wake  centerline.  On 
the  solid  surface  the  no-slip  condition  was  Imposed  along  with  the  second 
normal  derivative  of  pressure  set  to  zero.  At  the  outer  boundary  ellipse 
the  normal  second  derivatives  of  total  pressure  and  cartesian  transverse 
velocity  were  set  to  zero.  The  static  pressure  on  this  boundary  was  either 
set  to  the  value  determined  by  potential  flow  about  the  5:1  ellipse  (Zone 
Embedding  without  Interaction)  or  to  the  pressure  determined  by  the  inter¬ 
active  coupling  relation  to  be  discussed  below. 

Specification  of  Inflow  boundary  conditions  is  slightly  more  complex 
since  the  Inflow  properties  are  specified  through  the  boundary  layer  to  the 
wall.  The  boundary  conditions  are  the  specification  of  total  pressure, 
transverse  velocity  and  the  zero  first  streamwlse  derivative  of  static 
pressure.  These  conditions  in  cooperation  with  the  specified  static  pressure 
on  the  outflow  boundary  set  the  angle  of  attack,  the  mass  flux  through  the 
computational  zone  and  allow  the  static  pressure  at  the  Inflow  plane  to  be 
determined  by  the  rest  of  the  flow  field.  Allowing  the  Inflow  static  pressure 
to  adjust  to  downstream  flow  conditions  is  necessary  to  allow  transient 


24.  Douglas,  J.,  and  Gunn,  J.  E.:  A  General  Formulation  of  Alternating 
Direction  Methods.  Numerlsche  Math.,  Vol.  6,  1964,  pp.  428-453. 
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pressure  vsves  to  propsgate  upstream  and  out  of  the  computational  domain. 

The  specification  of  total  pressure  rather  than  streamvlse  velocity  maintains 
physically  reasonable  inflow  properties  as  a  pressure  transient  passes 
through  the  Inflow  boundary. 

The  difference  between  the  known  total  pressure  and  the  computed  static 
pressure  is  the  dynamic  head.  However  near  the  wall  this  quantity  approaches 
zero.  If  no  further  modeling  is  done,  an  Inconsistency  develops  near  the 
wall  as  the  static  and  total  pressures  approach  the  same  values.  In  fact, 
even  though  the  wall  point  on  the  inflow  boundary  need  not  be  included  in 
the  calculations  since  it  is  a  corner  point,  the  static  and  total  pressures 
are  equal  in  the  corner.  The  static  pressure  boundary  condition  is  used  to 
set  a  total  pressure  on  the  wall  which  is  consistent  with  the  interior  of 
the  computational  domain.  The  inflow  total  pressure  in  the  boundary  layer 
is  recomputed  each  time  step  to  be  consistent  with  the  Interior  of  the 
computational  domain.  The  Inflow  total  pressure  in  the  boundary  layer  is 
recomputed  each  time  step  to  be  consistent  with  the  free  stream  total 
pressure,  the  wall  pressure  and  the  input  dynamic  head  distribution  in  the 
boundary  layer  as: 

PT  ■  (pt##-  p,)f(y/8)  +  p# 

where  f(y/fi)  »  (-  )2  is  specified  on  input. 

“e 

The  interaction  between  the  Navier-Stokes  zone  and  the  Potential  Flow 
zone  starts  with  the  specification  of  a  streamline  near  the  body.  In  the 
results  presented  below  the  displacement  streamline  was  used.  This  stream¬ 
line  is  found  by  integrating  vertically  from  the  wall  at  a  series  of 
streamvlse  locations,  until  a  preset  mass  flux  has  been  captured.  In  order 
to  make  this  Navier-Stokes  determined  streamline  a  streamline  of  the  potential 
flow,  the  superposition  principle  for  incompressible  potential  flow  is  used. 

A  series  of  mass  sources  is  placed  on  the  streamline  whose  strengths  are 
determined  to  keep  the  potential  flow  about  the  body  plus  the  flow  from  the 
sources  from  crossing  the  streamline.  An  additional  source  was  placed  on 
the  forebody  to  account  for  the  flow  at  the  mid-chord  between  the  interaction 
streamline  and  the  ellipse.  The  net  effect  of  this  source  distribution  plus 
the  potential  flow  about  the  body  results  in  the  potential  flow  about  the 
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body  results  in  the  potential  flow  about  the  streamline .  This  new  potential 
flow  solution  determines  the  static  pressure  boundary  condition  used  on  the 
outflow  boundary  of  the  Navier-Stokes  zone.  The  interaction  procedure  was 
found  to  require  negligible  computer  time  and  was  repeated  after  each  time 
step  in  the  Navier-Stokes  solution  procedure. 


RESULTS 


Interactive  Zone  Embedding  was  used  to  analyze  the  flow  about  a  rounded 

trailing  edge.  Three  test  cases  were  run  for  laminar  flow  over  the  rear  of 

a  5:1  elliptic  cylinder  at  a  frees tr earn  Mach  number  of  0.2  and  a  Reynolds 

3 

number  based  on  chord  of  10  .  Two  computational  domains  were  used  for  the 
embedded  Nqvier-Stokes  calculations.  The  larger  domain  extended  approxi¬ 
mately  one  chord  from  the  body,  see  Fig.  5,  and  the  smaller  domain  extended 
only  about  one  half  chord.  With  the  larger  domain,  calculations  were  run 

both  with  and  without  interaction  to  verify  that  since  the  interaction 

2 

pressure  correction  falls  off  as  1/r  ,  there  was  only  a  very  small  influence 
of  the  boundary  layer  and  wake  on  the  static  pressure  at  the  outer  computa¬ 
tional  boundary. 

The  interacted  solutions  in  the  larger  and  smaller  domains  were  then 
compared.  An  overview  of  these  two  solutions  is  first  presented  by  comparing 
the  solution  in  the  smaller  domain  with  the  solution  in  the  corresponding 
portion  of  the  larger  domain  shown  in  Fig.  6.  Contours  of  constant  horizontal 
and  vertical  components  of  velocity,  vorticity,  and  total  and  static  pressure 
are  compared  in  Figs.  7-11.  Figure  12  shows  the  velocity  vector  comparison. 
Particularly  Insofar  as  the  recirculation  region  is  concerned  the  agreement 
between  the  larger  and  smaller  domain  calculations  is  very  good. 

As  expected,  the  static  pressure  is  the  most  sensitive  parameter  and, 
as  shown  in  Fig.  11,  it  exhibited  the  most  variation.  The  most  important 
part  of  the  static  pressure  field  is  the  pressure  distribution  at  the  airfoil 
surface,  since  this  can  be  integrated  to  calculate  force  and  moment  coeffi¬ 
cients.  The  wall  static  pressure  coefficients  for  the  larger  domain,  smaller 
domain,  and  for  the  reference  potential  flow  are  presented  in  Fig.  13.  The 
results  of  the  larger  and  smaller  domain  are  seen  to  be  virtually  identical. 

The  interaction  streamlines  for  the  two  domains  are  presented  in  Fig.  14. 
These  are  nearly  Identical  except  at  the  downstream  end  of  the  smaller  domain 
where  the  streamline  shape  was  extrapolated.  Outside  the  computational  domain 
the  interaction  streamline  was  assumed  to  follow  a  potential  flow  streamline. 
The  current  treatment  of  the  displacement  streamline  at  and  beyond  the  compu¬ 
tational  domain  is  suspected  to  be  one  source  of  the  discrepancy  between  the 
larger  and  smaller  domain  results  away  from  the  body.  Figure  15  shows  that 


the  velocity  and  static  pressure  on  the  wake  centerline  for  the  two  calcu¬ 
lations  are  in  good  agreement  near  the  body  but  that  this  agreement  deter¬ 
iorates  downstream  near  the  outflow  boundary. 

Traditionally  Interactive  Zone  Embedding  had  been  run  with  a  boundary 
l*y«t  analysis  for  the  inner  zone.  The  displacement  body  was  used  for  the 
interaction  since  it  is  both  a  streamline  of  the  flow  and  a  readily 
identifiable  boundary  layer  parameter.  However,  in  the  boundary  layer 
application  the  static  pressure  at  the  displacement  body  is  the  same  as  at 
the  boundary  layer  edge,  (normal  pressure  gradients  negligible  in  the 
boundary  layer)  but  this  is  not  necessarily  an  accurate  approximation  in 
the  present  case.  Since  the  shear  layer  is  not  thin  there  and  the  geometry 
nontrivial  there  are  significant  static  pressure  variations  beyond  the 
displacement  body,  and  the  integration  of  the  normal  pressure  from  the 
displacement  body  to  the  outer  computational  boundary  might  not  be  the 
same  for  the  (assumed)  potential  flow  as  it  is  for  the  Navler-Stokes  equations. 
This  could  contribute  to  the  static  pressure  at  the  outflow  boundary  of  the 
waller  computational  domain  falling  between  the  values  computed  by  the 
potential  flow  and  the  larger  domain  Navler-Stokes  at  the  same  locations  as 
shown  in  Fig.  16.  Note  that  Fig.  16  also  shows  a  kink  in  the  smaller  domain 
static  pressure  in  the  region  of  the  extrapolated  streamline  as  discussed 
in  the  previous  paragraph.  The  effects  of  choosing  a  streamline  much  closer 
to  the  outer  boundary  of  the  viscous  region  as  the  inner  boundary  for  the 
potential  flow  requires  further  investigation.  The  treatment  of  the  outflow 
boundary  Including  the  extrapolation  of  the  interaction  streamline  could 
also  be  refined. 


SUMMARY  AND  CONCLUSIONS 


The  Interactive  Zone  Embedding  concept  is  a  generalization  of  concepts 

which  have  been  used  for  some  time  with  boundary  layer  analyses.  This 

concept  has  been  extended  to  more  complex  flows  using  a  potential  flow  for 

the  outer  zone  and  the  compressible  Navier-Stokes  equations  for  the  inner 

zone.  In  this  interactive  manner  the  flow  about  a  rounded  trailing  edge 

was  computed  for  two  locations  of  the  outer  boundary.  The  calculations  were 

performed  on  a  5:1  elliptic  cylinder  at  a  Mach  number  of  0.2  and  a  Reynolds 

3 

number  based  on  chord  of  10  .  A  body  fitted  conformal  coordinate  system  was 
used.  In  the  first  Instance  the  outer  boundary  was  located  approximately 
one  chord  from  the  ellipse  and  in  the  second  the  outer  boundary  was  brought 
to  about  one-half  chord  of  the  ellipse. 

The  displacement  body  was  chosen  as  the  interactive  streamline  and  the 
potential  flow  over  this  streamline  was  computed  as  the  sum  of  the  analytic 
potential  flow  over  the  elliptic  cylinder  plus  the  flow  from  a  series  of 
sources.  The  strengths  of  these  sources  were  determined  to  make  the  inter¬ 
action  streamline  also  be  a  streamline  in  the  potential  flowfield.  The 
resulting  potential  flow  was  used  to  evaluate  the  static  pressure  at  the 
outflow  boundary  of  the  Navier-Stokes  zone. 

The  calculations  were  performed  with  the  large  Navier-Stokes  domain 
and  were  repeated  with  a  smaller  Navier-Stokes  domain.  The  flow  properties 
in  the  vicinity  of  the  body  for  the  two  calculations  were  in  excellent 
agreement.  Farther  from  the  body  the  fluid  properties  were  in  somewhat 
poorer  agreement.  It  is  anticipated  that  refinements  in  the  interaction 
model  could  Improve  the  flowfield  prediction  away  from  the  body.  In  partic¬ 
ular,  it  is  felt  that  locating  the  inner  boundary  for  the  potential  flow 
calculation  further  from  the  body  and  taking  care  not  to  bring  the  downstream 
boundary  of  the  viscous  region  too  close  to  the  ellipse  trailing  edge  would 
be  very  beneficial. 
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Figure  2.  Effect  of  variation  of  rear  stagnation  point  on  pressure 
distribution  predicted  according  to  Gostelov  (Ref.  3). 
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Figure  3.  Deviation  angle  comparisons  of  Ref.  4  reproduced 
from  Gostelow  (Ref.  3). 


A  analysis 


J  •'formation  provided  by  coupling  relations 


Figure  7.  Comparison  of  horizontal  velocity  contours 
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Figure  9.  Comparison  of  vortldty  contours 


Figure  10.  Comparison  of  stagnation  pressure  contours 
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Smaller  domain 
solution 


Figure  12.  Comparison  of  velocities 
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Figure  13.  Comparison  of  computed  wall  pressure  coefficients 


